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Abstract 

We have used the technique of spectro-astrometry to study the milli-arcsecond scale structure 
of the emission lines in the T Tauri star RU Lupi. The wings of the Ha emission are found to 
be displaced from the star towards the south-west (blue wing) and north-east (red wing) with 
angular scales of 20-30 milli-arcsecs. This structure is consistent with a bipolar outflow from 
the star. From a study of the variability of the intensity and position spectra we argue that 
a combination of magnetically-driven bipolar outflow and accreting gas contribute to the Ha 
emission. On the other hand, the [01] and [SII] emission are displaced from the star to the 
south-west but at much larger distances than the Ha, hundreds of milli-arcsecs for the high- 
velocity component (HVC) and down to 30 milli-arcsecs for the low- velocity components (LVC). 
The presence of both red-shifted and blue-shifted outflows in Ha but only a blue-shifted outflow 
in the forbidden lines can be explained if the disc obscures the red-shifted forbidden line outflow, 
but a disc gap with outer radius 3-4 au allows the red-shifted Ha to be seen. This gap could be 
induced by an unseen companion. 
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1 Introduction 



Optical emission lines from pre-main-sequence 
stars are observed as clues to understanding the 
acceleration of jets and mass accretion on to 
the stars. The spatial resolution of most tech- 
niques is insufficient to resolve their dynam- 
ics directly, while high-resolution spectroscopy 
has provided a variety of information for these 
mechanisms. Ha emission shows a rich variety 
of profiles due to an interplay between emis- 
sion and absorption features (e.g., Reipurth, 
Pedrosa, & Lago 1996), and various models 
have been discussed to explain their profiles 
using accelerating stellar wind (e.g., Decam- 
pli 1981), decelerating stellar wind (Mitskevich, 
Natta, & Grinin 1993), disc wind (Calvet, Hart- 
mann, & Hewett 1992), and infalling envelope 
(e.g., Calvet & Hartmann 1992). The forbid- 
den lines often have two blue-shifted velocity 
components, which have provided constraints 
on geometry and mechanism of jet and/or wind 
(Eisloffel et al. 2000, and references therein). 
The lack of red-shifted components of these line 
supports the existence of circumstellar discs (Eisloffel 



Lupi, is one of the most active T-Tauri stars in 
the southern sky. The Ha emission has broad 
wings up to 900 km s _1 (Reipurth, Pedrosa, 
k, Lago 1996), and the equivalent width of the 
emission line is up to 200 A, which is one of 
the largest among T-Tauri star (Giovannelli et 
al. 1995). Forbidden lines of [SII] and [OI] are 
also present in the spectrum of RU Lupi, and 
at least the [OI] lines have an extended blue 
wing up to 300 km s _1 and a redder peak up 
to 30 km s _1 (Gahm, Lago, & Penston 1981; 
Hamann Sz Persson 1992). In this paper, the 
detail of the observations and the data reduc- 
tion are described in §2, the results are shown 
in §3, and the origin of Ha emission and the for- 
bidden lines are discussed in §4 together with 
the existence of a circumstellar disc with a gap. 
Throughout this paper, the distance of 140 pc 
to the Lupus cloud is adopted based on Hughes, 
Hartigan, & Clampitt (1993). 



Observations & Data Reduc- 
tion 



et al. 2000, and references therein). In addi- 
tion to high-resolution spectroscopy, some tech- 
niques have been tried in the last decade to ob- 
tain spatial information. Speckle observations 
reveal elongated structure of the Ha emission 
with an angular scale of 0.09 arcsec towards T 
Tau (Devaney et al. 1995). The Hubble Space 
Telescope has provided fine structure of jets to- 
wards young stellar objects (Ray et al. 1996; 
Burrows et al. 1996) although the bandwidths 
of its filters are often insufficient to suppress 
the stellar continuum from pre-main-sequence 
stars. Ground-based long-slit spectroscopy is 
effective in measuring the typical angular scale 
of a few forbidden lines suppressing the con- 
tamination of the strong continuum emission 
(e.g., Hirth, Mundt, and Solf 1997). 



Observations were carried out on 1996 Aug 25, 
1997 Jun 27, and 1999 Jul 2 at the 3.9-m Anglo- 
Australian Telescope using the RGO spectro- 
graph. The configuration with a 1-arcsec slit 
width, a 1200 line mm -1 grating, and the 82- 
cm camera provided a spectral resolution (A/AA) 
of 7000. The first two observations used a Tek- 
tronix 1024x1024 thinned CCD. The 1999 Jul 
2 observation used a 2048 x 4096 deep depletion 
CCD from MIT Lincoln Labs. The pixel scale 
was 0.23 arcsec with the Tektronix CCD and 
0.15 arcsec with the MITLL CCD, which pro- 
vides good sampling of the seeing profile (1-1.5 
arcsec) at each wavelength. The spectra were 
obtained at four slit position angles (0°, 90°, 
180°, and 270°). The flat fields were made by 
combining many exposures with the spectro- 



Here we present the results of spectro-astrometiSgaph illuminated by a tungsten lamp. Wave- 
observations of RU Lupi. The spectro-astrometry length calibrations were made by observations 
technique studies the small scale structure by °f a CuAr lamp. 

observing the relative position of the centroid The data were reduced using the FIGARO 
of the point-spread function on scales of a few package. After subtracting the bias level and 
milli-arcsec (Bailey 1998a, b). The target, RU flat-fielding, the position spectrum was deter- 
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mined by fitting the seeing profile at each wave- 
length with a Gaussian function. About 1- 
2xl0 4 photons at each wavelength provides a 
typical accuracy of 2-5 milli-arcsec at the con- 
tinuum level. Any instrumental effect in the 
position spectra were eliminated with high ac- 
curacy by subtracting those with opposite po- 
sition angles (0 ° - 180°, or 90 ° - 270°), since 
these will remain constant on rotation whereas 
any real structure in the source will reverse its 
sign. Effects which are removed in this way in- 
clude any curvature or optical distortion intro- 
duced by the spectrograph, any misalignment 
of the spectrum with the CCD columns, and 
any effects due to departure of the CCD pix- 
els from a regular grid, imperfect flat fielding, 
or charge transfer deficiencies in the CCD. The 
derived position spectra have wavelength cov- 
erages of 6520-6750 A for the former two runs 
and 6160-6780 A for the latest run, which in- 
clude Ha, [SII] 6716/6731 A, [Nil] 6583 A, and 
[OI] 6300 A. Each position spectrum has an 
arbitrary zero point, which we adjust to cor- 
respond to the continuum position. In addi- 
tion to the position spectra, the intensity spec- 
tra were obtained by subtracting the bias, flat- 
fielding, subtracting the adjacent sky and ex- 
tracting bright columns on the CCD. In all the 
spectra, the radial motion of the star is cali- 
brated using the Li 6707. 815A absorption line 
on the stellar atmosphere. The measured ra- 
dial velocity in the local standard rest frame is 
+8 + 2 km s _1 . 



sity spectra, a number of emission lines includ- 
ing Ha, Hel, [SII], [OI], [Nil], Fel, Fell, and 
Sill are detected. The wavelengths and equiva- 
lent widths without atmospheric correction are 
summarized in Table 1. 

Among these lines, the Ha, [SII] 6731/6716A, 
[Nil] 6583 A, and [OI] 6300 A lines show dis- 
placements from the continuum position in the 
position spectra. The positions shown in Fig. 1 
are the centroid of the combined line and con- 
tinuum emission of the star, that is, the cen- 
troid of their positions weighted by their inten- 
sities. To get the true position of the emission 
line region the observed position must be cor- 
rected for the effect of the underlying contin- 
uum by dividing the positional displacement in 
the raw spectra by I A (iinc)/(^A(iinc) + -^A(conn) ) ■ 
Fig. 2 shows the positional displacement of 
Ha, [OI] 6300 A, and [SII] 6731/6716 A ob- 
served on 1999 Jul 2 after this calibration. Typ- 
ical uncertainties for the displacements in Fig. 
2 are 2, 15, 40, and 100 milli-arcseconds for 
Ha, [OI] 6300 A, [SII] 6731 A and 6716 A, re- 
spectively. For Ha, higher accuracy was ob- 
tained than that at the continuum level be- 
cause of its brightness while the accuracy of 
the other forbidden lines are much less than 
the continuum because of their lower intensity. 
The line profiles shown in Fig. 2 are similar to 
those previously presented by Giovannelli et al. 
(1995), Reipurth, Pedrosa, & Lago (1996), and 
Hamann (1994). 

Fig. 2 shows that the positional displace- 



Similar methods were previously used by Hirth,ment of Ha emission extends towards the south- 



Mundt, and Solf (1994, 1997) to study forbid- 
den emission lines in pre-main sequence stars. 
However, the techniques described here enable 
us to achieve substantially higher accuracy in 
the relative positions. 

3 Results 

3.1 Positional displacement and in- 
tensity profile of Ha, [OI] and [SII] 
emission 

Fig. 1 shows the raw spectra of the intensity 
and the positional displacement. In the inten- 



west and the north-east over angular scales of 
up to 20-30 milli-arcsec. The intensity pro- 
file of Ha emission is symmetric except for a 
redshifted absorption dip at +50 km s~ 1 from 
the star, and the blue-shifted wing corresponds 
to the displacement at the south-west while 
the red-shifted wing corresponds to that at the 
north-east. At these wings, the displacement 
monotonically increases with the relative veloc- 
ity from the star in the velocity ranges of —50 
to -200 km s" 1 and +200 to +300 km s -1 . On 
the other hand, the displacement is distributed 
within 5 milli-arcsec from the star in a lower 
velocity range between —50 and +200 km s _1 , 
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and does not show clear correlation with the 
velocity in Fig. 2. 

The positional displacement of the three for- 
bidden lines extends only towards the south- 
west, the same direction as those of the blue- 
shifted Ha wing, and that of HH55 which lies 3 
arcmin from RU Lupi (cf. Krautter, Reipurth, 
& Eichendorf 1984). The displacements arc 
larger than that of Ha emission by a factor of 
10-30, and smallest in [01] and largest in [SII] 
6716A: up to 200, 400, and 600 milli-arcsecs 
for [OI] 6300A, [SII] 6731 A, and 6716A, re- 
spectively. These line profiles show two peaks 
with typical velocities between —250 and —100 
km s _1 for the high velocity component (HVC) 
between —100 and km s _1 for the low velocity 
component (LVC). The HVC is more extended 
spatially than in the LVC as in other T Tauri 
stars. The typical angular scales of the LVC are 
30-40 milli-arcsec in [OI] 6300 A and 150 milli- 
arcsecs in [SII] 6731 A. Fig. 2 does not show 
a clear difference in the displacement between 
the two components in [SII] 6716A. 

The relation between the velocity and the ra- 
dial distance of the displacement for Ha, [OI] 
6300A, and [SII] 673lA, is plotted in Fig. 3. A 
higher signal to noise ratio of the positional dis- 
placements has been obtained by binning along 
the wavelength. This figure shows that the 
velocity monotonically increases with the dis- 
tance for all the lines, and the gradient for Ha 
emission is much larger than those for the for- 
bidden lines. Fig. 3 also shows slight asym- 
metry of the velocity field between the blue- 
shifted and red-shifted wings of Ha emission: 
the positional displacement is prominent in the 
figure beyond the velocity of 100 km s -1 for the 
blue-shifted wing and 200 km s _1 for the red- 
shifted wing. In the forbidden lines, the dis- 
played velocity gradient is larger in the HVC 
(AV > 100 km s _1 ) and smaller in the LVC 
(AV<100 km s -1 ). On the other hand, the 
LVC in [OI] reveals no velocity gradient while 
those in the [SII] 6731 A line show slight posi- 
tive gradients. The position of the LVC is offset 
from the star even at the zero velocity: 30 and 
70 milli-arcsecs for [OI] and [SII] 6731 A, re- 
spectively. 



The observed tendencies for the line profiles 
and the positional displacement of the HVC 
agree with the previous observations of other 
T Tauri stars. Ground-based long slit spec- 
troscopy shows the positional displacement of 
the HVCs of T Tauri stars in the Taurus- Auriga 
cloud (Hirth, Mundt, and Solf 1997), which lies 
at a similar distance (140 pc) to that of the Lu- 
pus cloud. The HVC and LVC in these objects 
have typical velocities —50 to —150 km s^ 1 and 
—5 to —20 km s _1 , respectively, which are sim- 
ilar to those of RU Lupi. The measured dis- 
placements of the HVC are 0.2 and 0.6 arcsecs 
for [OI] and [SII] 673lA, respectively, which 
are similar to those seen in our observations. 
The [SII] lines have much lower critical density 
than the [OI] line (typically 10 3 and 10 6 cm -3 , 
respectively), and the observed positional dif- 
ference between these lines indicates that the 
outer region has a lower electron density. Such 
a distribution of the electron density could be 
due to diverging stream lines of the outflow, as 
suggested by Hirth, Mundt, and Solf (1997). 

3.2 Time variation of Ha and [SII] 
emission 

Fig. 4 shows the time variation of the intensity 
profile and the positional displacement for Ha, 
[SII] 6731A, and [SII] 6716A emission. Since 
photometric observations were not made on the 
same dates, each line profile is displayed in a 
unit of the continuum intensity on each date. 
The positional displacements are displayed af- 
ter removal of the contamination of the contin- 
uum, as in Fig. 2. 

In the case of Ha, time variation was de- 
tected for both the line profile and the posi- 
tion spectra. The intensity profiles observed 
on 1996 Aug 25 and 1997 Jun 27 are nearly 
the same at the wings. On the other hand, at 
the peak with the velocity of —100 to +100 km 
s _1 , the profile on the latter date has lower in- 
tensity relative to the continuum than that on 
the former date. On 1999 Jul 2, the relative 
intensity of the Ha to the continuum is slightly 
lower than the others except that at the ab- 
sorption dip. In the same observing run, the 
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absorption dip is shallower and relative inten- 
sity is higher at the top than those observed 
in the other runs. Clear change of the position 
spectra of Ha was observed during 1997-1999 
with the velocity ranges of —200 to —100 km 
s _1 and +100 to +300 km s _1 . The positional 
displacements of these wings are larger on the 
latest date than in the others by a factor of 
2-3. The repeatability of the measured posi- 
tional displacement between —100 and +100 
km s _1 confirms that the displacement is not 
due to slight oscillation of the tracking of the 
telescope, since any other pre-main sequence 
stars observed in 1999 do not show the same 
displacements as observed in this object. 

The relative intensity of the [SII] lines to the 
continuum is smallest on the earliest date, and 
largest on the latest date at the whole velocity 
range. On the other hand, the peak velocity 
and the line widths of the HVC and LVC did 
not change throughout the observations. To 
determine the difference of the time variation 
between the HVC and LVC, we measure the 
equivalent widths of the HVC and LVC respec- 
tively, and these are tabulated in Table 2. For 
both the [SII] 6716A and 673lA, the LVC shows 
a larger increase in the equivalent width than 
the HVC: the difference of the relative inten- 
sities between 1996 and 1999 is 40-50 and 10- 
30 percent for the LVC and HVC, respectively. 
Table 2 also shows the equivalent width ratio of 
the two lines of 0.53-0.60 and 0.37-0.40 for the 
HVC and LVC, respectively, and these ratios 
indicate electron densities of 4-6 x 10 3 and more 
than 10 4 -10 5 cm 3 assuming a electron temper- 
ature of 10 4 K (cf. Osterbrock 1989). Such a 
difference of the electron density between the 
HVC and LVC was also observed in other T 
Tauri stars by Hamann (1994). On the other 
hand, Fig 4. shows no systematic time vari- 
ation in the position spectra of the two [SII] 
lines. The displayed variation would indicate 
motion of gas, although higher signal-to-noise 
ratio is required to comfirm the existence of the 
variation. 

In addition to the time variation of the emis- 
sion line fluxes, variation of the continuum flux 
can affect the intensity profiles shown in Fig. 



4. Its contribution is not clear owing to the 
lack of photometric data. Previous observa- 
tions show the time variation of the contin- 
uum which may be caused by obscuration by 
clumpy clouds (Gahm et al. 1974; Giovannelli 
et al. 1995). In addition, Table 1 shows time 
variation of the equivalent width for the Li ab- 
sorption, which suggests variation of the veiling 
continuum. The intensity profiles of Ha and 
the two [SII] lines show different time variation, 
and it cannot be simply explained by variation 
of the continuum. The ratios of the equiva- 
lent width between the LVC and the Li absorp- 
tion are nearly constant over time (1.4-1.6 and 
0.5-0.6 for [SII]673lA/Li and [SII]6716A/Li, re- 
spectively), and suggest that the time variation 
of the LVC is due to the variation of the veiling 
continuum. 

3.3 Comparison with the other lines 
profiles 

The profiles of the other lines are plotted in 
Fig. 5. Most of the permitted lines includ- 
ing Fel, Fell, Sill have triangular profiles with 
peak velocities of —50 to +50 km s _1 , and full 
width zero intensities of less than 250 km s -1 . 
The width of these lines is much narrower than 
that of the Ha, and the velocity ranges of these 
lines do not overlap that of the Ha wings: thus, 
these lines are not related to the Ha wings ob- 
served in 1999. On the other hand, the profiles 
of [OI] 6364 A and [Nil] 6583 A are much more 
similar to those of [OI] 6300 A and the two [SII] 
lines than that of Ha. The profiles of [OI] 6364 
A, which shares the upper transition level with 
[OI] 6300 A, is nearly the same as that of the 
latter line shown in Fig. 2. No clear signature 
of the LVC is found in the [Nil] profile as in 
the case of other T Tauri stars (e.g., Hartigan, 
Edwards, and Ghandour 1995). 

The Hel 6678 A line shows different profiles 
from the others, and a prominent time varia- 
tion was observed. In 1996 and 1997, the line 
had nearly the same triangular profile with a 
full width zero intensity of about 400 km s _1 . 
On the other hand, in 1999, the line width 
was reduced to be about a half without chang- 
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ing the peak velocity and the relative inten- 
sity to the continuum at the peak. Two com- 
ponents with different line widths exist in the 
latter profile: a narrow component (NC) with 
a full width half maximum of about 50 km 
s _1 , and a broad component (BC) with a full 
width zero intensity of more than 200 km s _1 . 
Such coexistence of two components has been 
often observed in Hel 5876A, Call, and Fell 
lines of other T Tauri stars (e.g., Batalha et 
al. 1996, Muzerolle, Hartmann, & Calvet 1998, 
Beristain, Edwards, & Kwan 1998), and both 
the components are considered to be associ- 
ated with mass accretion in the vicinity of the 
star (e.g., Najita et al. 2000, and references 
therein). As shown in Fig. 4, Ha profiles have 
a red absorption dip indicative of infalling ma- 
terial to the star. The Ha and Hel profiles ob- 
serving in 1999 are characterized by a shallower 
redshifted absorption dip and a weaker broad 
component, respectively, and such relation be- 
tween the Ha dip and the Hel broad component 
suggests that the latter component arises from 
gas accreting to the star. No prominent emis- 
sion associated with the Ha wings was found 
in the broad component of Hel in 1999: thus, 
we conclude that the Ha wings have a different 
origin from that of the broad component of the 
Hel line. 



4 Discussion 

4.1 Origin of Ha emission 

Positional displacement as measured by spectro- 
astrometry should be sensitive to any asym- 
metric structure about the continuum source. 



flow. No binary companion has been detected 
by the Hubble Space Telescope (Bernacca et 
al. 1995) or infrared speckle observations (Ghez 
et al. 1997). Thus, it is likely that the detected 
feature reflects the distribution of circumstellar 
matter around the star. 

A variety of mechanisms have been previ- 
ously proposed to explain the observed Ha pro- 
files from pre-main-sequence stars. Magnetically- 
driven stellar winds which accelerate with dis- 
tance have been suggested for Balmer emission 
including Ha (e.g., Decampli 1981; Hartmann, 
Edwards, & Avrett 1982; Lago 1984; Natta, 
Giovanardi, &i Palla 1988), while a decelerat- 
ing "stochastic" stellar wind was proposed by 
Mitskevich, Natta, & Grinin (1993). Calvet 
& Hartmann (1992) and Hartmann, Hewett & 
Calvet (1994) argue that the Ha emission is 
produced in infalling envelopes. Edwards et al. 
(1994) presented for 15 T Tauri stars a set of 
high resolution profiles including Balmer, Hel, 
and Nal lines and showed that the blue/red 
asymmetry of upper Balmer lines agrees well 
with magnetospheric accretion models. On the 
other hand, Reipurth, Pedrosa, and Lago (1996) 
compared their Ha profiles with various mod- 
els, and suggested that a single model cannot 
account for the variety of the Ha profiles. In 
addition to outflowing and inflowing motion, 
Keplerian rotation would also be responsible 
for the shape of the Ha profile. 

Among these candidates, we consider bipolar 
outflow as a plausible explanation for the ob- 
served displacement for the following reasons: 
(1) the positional displacement of the blue-shifted 
component is in the same direction as those of 
forbidden lines, which are considered to trace 
In the case of Ha emission from pre-main-sequencejutflow, and (2) the other permitted lines, which 



stars, binary companions are responsible for 
the displacement as shown by Bailey (1998a). 
However, the feature seen in RU Lupi is quite 
unlike those normally seen in binary systems 
which are not symmetric about the line centre. 
Binary stars with the same spectrum and differ- 
ent radial velocities could produce the observed 
positional displacement, although such motion 
at the observed direction does not agree with 
the orbital motion perpendicular to the out- 



probe mass accretion, do not have such broad 
wings as Ha, thus these wings should have dif- 
ferent origin (cf. §3-3). No clear evidence for 
inflowing motion has been detected in the po- 
sition spectra, consistent with previous models 
which propose that the Ha emission arises from 
a region within less than ten stellar radii (e.g., 
Calvet & Hartmann 1992; Hartmann, Hewett 
& Calvet 1994). Rotation cannot explain the 
observed displacement for the following reasons: 
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(1) the measured velocity increases with po- 
sitional displacement while Keplerian rotation 
should behave in an opposite fashion, and (2) 
the star needs to be more than 100 M Q to 
bound the gas with the measured positional 
displacement and velocity (3 au and 200 km 
s _1 , respectively) while previous research sug- 
gests the mass of the star is 0.5-2 M (Gahm 
et al. 1974; Lamzin et al. 1996). 

Figs 2 and 3 show that the velocity of the 
outflowing gas increases with positional displace- 
ment, a tendency which agrees with that of 
magnetically-driven wind models which predict 
that the flow accelerates with distance (e.g., 
Decampli 1981; Hartmann, Edwards, & Avrett 
1982; Lago 1984). Most of the models adopt a 
spherically symmetric line forming region not 
exceeding ten stellar radii, and predict that the 
acceleration of the flow ends within this scale. 
On the other hand, the observed position spec- 
tra of RU Lupi indicate bipolar geometry and a 
velocity gradient extending to a few au (Figs 2 
and 3). The latter results would indicate slower 
acceleration of the outflow than those predicted 
by models, however, models with more appro- 
priate geometry and direct prediction of the po- 
sitional displacement are necessary for detailed 
comparison. Actually, Ha emission from pre- 
main-sequence stars should be optically thick 
(cf. Reipurth, Pedrosa, Sz Lago 1996), thus the 
observed velocity field would be highly sensitive 
to the geometry which affects radiative trans- 
fer. 

The time variation of the positional displace- 
ment and the intensity profile for the Ha wings 
are similar to that of Hel in the following points: 
there is no clear difference between 1996 and 
1997, and a remarkable difference between 1997 
and 1999. Since the Hel emission is related to 
the accreting g shown in §3-3, such similar 
time variation suggests that the accreting gas 
also contributes to the Ha flux at least in 1996 
and 1997. Even if the positional displacement 
of the bipolar outflow is constant over time, the 
net positional displacement can change as the 
flux of the accreting gas at the centre varies. 
Transient time variation of the displacement in 
high-velocity wings could be explained by di- 



rect motion of the outflowing ionized gas, al- 
though we reject this possibility since line pro- 
files with the same velocities only show slight 
change of line-to-continuum ratio at the same 
velocities. 

Here we attempt to explain the observed po- 
sitional displacement with a steady bipolar out- 
flow and time- variable flux from inflow. To re- 
move contamination from inflow on the total 
Ha flux, and reproduce the position spectra 
of the outflow we make the the following as- 
sumptions: (1) the Ha emission from infalling 
gas arises from the star and cannot be spa- 
tially resolved by our observations, (2) inflow 
contributes to the observed Ha profiles only 
in 1996 and 1997, and (3) the time variation 
of the Li equivalent widths shown in Table 1 
is due to the veiling continuum. The first as- 
sumption provides the observed net positional 
displacement as follows: 

Ax x = Ax Xout ■ - /A(out) (1) 

I X(out) + i A(in) 

where Ax\ is the net displacement shown in 
Fig. 4, AxXout is the positional displacement 
caused by the bipolar outflow, iu ou t) and ^A(in) j 
are the absolute Ha intensity from the out- 
flow and inflow, respectively. This assumption 
is consistent with our observations, and the 
previous inflow models described above. From 
the second assumption, the intensities from the 
outflow and inflow are obtained as follows: 

^A(out) = ^A(Hq99)j I\(in) = ^A(Ha) ~ A(Ha99)> 

(2) 

where I\m a ) 1S the absolute Ha intensity in a 
certain observing run, and /A(Ha99) is that ob- 
served in 1999. This assumption is plausible for 
the Ha wings since those observed in 1999 are 
not associated with the other permitted lines 
indicative of mass accretion as described in 3- 
3. The third assumption provides the absolute 
Ha intensities from the observed parameters as 
follows: 

h(Ha) i\(Ha) [ -^A(conn) _ ^A(Ha) ^A(Li99) 

^A(Ho;99) h(Ra99) 1 ^A(conn99) *A(Ha) W X (Li) 

(3) 
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where i\m a ) an d ^A(conn) are the relative Ha 
intensity to the continuum and the absolute in- 
tensity of the continuum in a certain observing 
run, respectively, and i A (Ha99) and ^A(conn99) are 
those observed in 1999. Wxa,i) and HA(Li99) are 
the equivalent widths of the Li absorption in a 
certain observing run and in 1999, respectively. 
By combining equations (l)-(3), the positional 
displacement of the outflow is obtained from 
the observed parameters as follows: 

a _ A ™ ^(Hn) W X(Li99) 
As A(out) - Asa • 777 , (4) 

*A(Hq99) " A(Li) 

Fig. 6 shows the reproduced position spec- 
tra using equation (4). The spectra for 1999 
are identical with the observed one shown in 
Fig. 4, since we assume that the outflow dom- 
inates the observed intensity in this observing 
run. These spectra agree well with each other, 
thus support the idea that the time variation 
of the positional displacement can be explained 
by the combination of steady bipolar outflow 
and time- variable flux from the accreting gas. 

The contribution of the inflow to the total 
Ha flux could be estimated from equations (2) 
and (3): according to these equations, 32 and 
28 percents of the total Ha flux would arise 
from the accreting gas in 1996 and 1997, re- 
spectively. These values would be lower lim- 
its, since the inflow with low velocities would 
contribute to the Ha emission in 1999. The 
permitted lines related to mass accretion were 
detected with a typical velocity range between 
— 150 and +150 km s _1 as shown in §3-3, and 
there is no reason to assume that Ha emission 
does not arises from the same region. 

4.2 Origin of the forbidden lines 

Optical forbidden emission lines including [01] 
6300A and [SII] 6716/6731 A provide informa- 
tion on the outflows from T Tauri stars. These 
line profiles often have two blue-shifted velocity 
components: a high-velocity component (HVC) 
and low-velocity component (LVC) which have 
typical velocities of —50 to —150 km s _1 and 
—5 to —20 km s -1 , respectively (Hirth, Mundt, 
and Solf 1997). The HVC is spatially resolved 



by ground-based long-slit spectroscopy in the 
nearest star- forming regions (e.g., Solf & Bohm 
1993; Hirth, Mundt, and Solf 1997). Hirth, 
Mundt, and Solf (1997) have shown that the 
typical angular scale of the HVC is 0.2-0.6 arc- 
sec from their large sample of data mainly in 
the Taurus-Auriga cloud. The high speed and 
narrow linewidth of the HVC indicates that the 
velocity vectors of the contributing gas parti- 
cles are nearly parallel to one another, i.e., sug- 
gesting a well-collimated flow or jet (Kwan & 
Tademaru 1988). On the other hand, the low- 
velocity component (LVC) has smaller angu- 
lar scales in the same sample (0.1-0.2 arcsec or 
less), and its origin is less understood. Geomet- 
ric effects of a single outflow are proposed to 
explain the existence of the two velocity com- 
ponents (e.g., Edwards et al. 1987; Hartmann 
&: Raymond 1989), while two distinct mech- 
anisms, jet and disc wind, are proposed for 
each velocity component by Kwan and Tade- 
maru (1988, 1995). Hamann (1994) measured 
the electron density, temperature, and the ion- 
ization state between the HVC and LVC in 
many T Tauri stars, and the measured differ- 
ence of the physical conditions suggest distinct 
origins for these two components as proposed 
by Kwan and Tademaru (1988, 1995). Sim- 
ilar results were also presented by Hartigan, 
Edwards, &; Ghandour (1995). The velocity 
dispersion and flow velocity are comparable in 
the LVC, suggesting moderate collimation such 
as that due to a disc wind (e.g., Solf Sz Bohm 
1993). Most pre-main-sequence stars do not 
have red-shifted components to their forbidden 
lines, and this is explained by obscuration of 
these components by circumstellar discs (e.g., 
Eisloffel et al. 2000, and references therein). 

The [OI] and [SII] lines from RU Lupi have 
two blue-shifted components, and our results 
also support the existence of two distinct mech- 
anisms. The [SII] intensity ratio indicates dif- 
ferent electron densities between the HVC and 
LVC as shown in §3-2 (4-6 xlO 3 and more than 
10 4 -10 5 cm~ 3 assuming a electron temperature 
of 10 4 K), a tendency also observed in other T 
Tauri stars by Hamann (1994). In addition, the 
intensity of each component exhibits different 
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time variations as described in §3-2. The coex- 
istence of two distinct mechanisms can easily 
explain these differences. 

Furthermore, the velocity field in the LVC 
shown in Fig. 3 agrees with the disc wind hy- 
pothesis. Fig. 3 shows that the position of 
the LVC is offset from the star even at zero ve- 
locity and nearly constant over all velocities in 
[01] emission, as described in §3-2. The former 
cannot be explained by acceleration of outflow, 
while both tendencies can easily be explained 
by moderate collimation and/or rotation which 
would be expected in a disc wind. In the case of 
[SII] 673lA, the positional displacement in the 
LVC increases with the velocity. Such discrep- 
ancy between the [01] and [SII] emission can 
be explained by blending of the HVC and LVC 
near their boundary. The profiles of the forbid- 
den lines in Fig. 2 show that the two compo- 
nents are partially blended, thus contamination 
from the HVC could increase the positional dis- 
placement in the LVC near their boundary. It 
is likely that this effect is more prominent for 
the [SII] line, since the lower critical density of 
this line allows for larger contamination from 
the HVC which has a lower electron density 
(cf. §3-2). 

4.3 Circumstellar disc and disc gap 

As shown in Fig. 2, the forbidden lines have 
only blue-shifted components, similar to many 
other T Tauri stars, a fact usually explained by 
obscuration of the red-shifted flow by a circum- 
stellar disc (e.g., Eisloffel et al. 2000, and ref- 
erences therein). Infrared excess and a bright 
millimetre continuum suggest the existence of 
a circumstellar disc around RU Lupi (e.g., Gio- 
vannelli et al. 1995; Carballo et al. 1992; 
Niirnberger et al. 1997), supporting the expla- 
nation for the lack of red-shifted components. 

On the other hand, Ha emission has both 
blue- and red-shifted components. The mea- 
sured positional displacement for the Ha emis- 
sion is smaller than any velocity component of 
[OI] and [SII] emission. The detection of the 
red-shifted Ha outflow can be explained by a 
gap in the disc which allows this emission to 



be seen. The presence of disc holes and gaps is 
often inferred around other pre-main-sequence 
stars by infrared spectral energy distributions 
(SEDs) (e.g., Mathieu et al. 1991; Marsh & 
Mahoney 1992; Jensen & Mathieu 1997), and 
directly observed in near-infrared imaging (e.g., 
Roddier et al. 1996; Silber et al. 2000) and at 
radio wavelengths (e.g., Dutrey, Guilloteau, & 
Simon 1994). Fig. 3 suggests that the pro- 
jected radius of the hole or gap is 20-30 milli- 
arcsec, which corresponds to 3-4 au, assuming 
symmetric distribution of the blue-shifted and 
unseen red-shifted components of the forbidden 
lines. 

To investigate the presence of a disc gap, we 
plot the infrared SED in Fig. 7. The flux 
at JHKLM-hands and 12, 25, 60, 100 u.m 
were obtained from Giovannelli et al. (1995) 
and Carballo et al. (1992). Giovannelli et 
al. (1995) show the near-infrared flux is highly 
variable with time, thus the data-set with the 
largest and smallest flux were selected from their 
data and plotted in the figure. In Fig. 7, 
a spectrum obtained via the Infrared Space 
Observatory (ISO) data archive is also plot- 
ted to interpolate the flux in the mid-infrared. 
Lamzin et al. (1996) claim that the mid-to-far 
infrared flux obtained by Carballo et al. (1992) 
would be overestimated, since another source 
discovered by Cohen & Schwartz (1987) is in- 
cluded in the same aperture. However, the con- 
sistency of their data with the ISO spectrum 
shows that the contribution from this source is 
much less than the total flux: the aperture of 
the ISO Short- Wave Spectrometer (SWS) with 
a size of 20x33 arcsec 2 does not include this 
source which lies about 3 arcmin away from 
the target. 

Fig. 7 shows two peaks in the SED which 
arises from dust. The near-infrared peak at 
1.65 /im is longer than the peak for the star 
(1.0 /Ltm) with an effective temperature of 3900 
K (Niirnberger et al. 1997; Lamzin et al. 1996). 
Previous authors show that the selective ex- 
tinction towards RU Lupi is 0.3-1.3 Ay (Hughs 
et al. 1994; Giovannelli et al. 1995). Such 
values can shift the peak of the stellar con- 
tinuum by less than 0.2 /im suggesting that 
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the near-infrared excess from the stellar pho- 
tosphere represents hot dust from the accre- 
tion disk. Such a contribution by the accretion 
disk to the near-infrared SED is also supported 
by the near-infrared colours (J-H and H-K of 
0.84-0.92 and 0.72-0.77, respectively) and dia- 
grams shown by Hughes et al. (1994), Greene 
& Meyer (1995), and Itoh et al. (1999). On the 
other hand, the mid-to-far infrared SED shows 
another excess which is due to cooler dust with 
a temperature of less than 100 K. The shallow 
dip at 4-15 /im suggests a lack of dusty material 
at temperatures of 200-900 K. These tempera- 
tures correspond to the radiative temperature 
at a radius of 0.1-2 au from the star, assuming 
the total luminosity from the star and accret- 
ing gas of 5 L Q based on Lamzin et al. (1996). 
These spatial scales agree well with that of the 
suggested gap in Figs 2 and 3. 

Theoretical work has shown that disc gaps 
and holes could be induced by binary compan- 
ions or young planets (e.g., Lubow & Artimow- 
icz 2000). In the case of RU Lupi, no other 
pre-main-sequence companions have been de- 
tected by Hubble Space Telescope (Bernacca et 
al. 1995), infrared speckle observations (Ghez 
et al. 1997), or spectro-astrometry which is 
sensitive to the presence of binaries (cf. Bai- 
ley 1998a). Therefore, the observed disc gap 
suggests the existence of an unseen compan- 
ion. A young planet could be one of possible 
candidates for the following reasons: the pres- 
ence of the fragmented dusty clouds around RU 
Lupi are suggested by variability at optical to 
infrared wavelengths (Gahm et al. 1974; Gio- 
vannelli et al.1995), and Gahm et al. (1974) 
claim that such clouds would indicate ongoing 
planetary formation around this object. In ad- 
dition, the observed disc gap has a similar size 
to that of the orbit of the Jovian planets in our 
solar system, and Takeuchi, Miyama, & Lin 
(1996) suggest that a proto-Jupiter could in- 
duce a disc gap with a similar size to its orbit. 
Their numerical calculations show that a Sat- 
urnian and Jovian mass protoplanet can open 
a gap in 10 3 -10 4 and 10 2 -10 3 yrs, respectively, 
timescales which are much shorter than the age 
of RU Lupi (10 6 yrs - Lamzin et al. 1996). 



5 Conclusions 

Spectro-astrometric observations have been used 
to determine the spatial distribution of optical 
emission lines from RU Lupi. Positional dis- 
placements were detected in Ha, [OI] 6300 A, 
[SII] 6731/6716 A, and [Nil] 6583 A with ac- 
curacies down to a few milli-arcsec. The inten- 
sity profiles were also obtained for these lines 
together with Hel, Fel, Fell, Sill, and another 
[OI] line. 

The positional displacement of Ha emission 
extends towards the south-west and north-east 
with angular scales of 20-30 milli-arcsecs. The 
blue-shifted wing in the intensity profile corre- 
sponds to the displacement at the south-west 
while the red-shifted wing corresponds to that 
at the north-east. We conclude that the posi- 
tional displacement of the Ha emission is due 
to bipolar outflow for the following reasons: (1) 
the distribution of the positional displacement 
is almost symmetric between the blue-shifted 
and red-shifted components, (2) the displace- 
ment of the blue wing aligns well with that 
of the forbidden lines which are considered to 
trace outflowing gas. The velocity of the out- 
flow increases with the distance, suggesting the 
flow is magnetically-driven. On the other hand, 
the time variation of the intensity profile and 
the positional displacement can be explained 
well by the variation of Ha flux from the ac- 
creting gas at the star. We estimate the con- 
tribution of the accreting gas to the total flux 
to be more than 32 and 28 percent in 1996 and 
1997, respectively. 

The positional displacement of [OI] and [SII] 
emission extends towards the south-west, in the 
same direction as HH55 which lies 3' from the 
star. The intensity profiles have two blue-shifted 
components, as seen in many other pre-main- 
sequence stars, and the positional displacements 
of hundreds of milli-arcsecs and down to 30 
milli-arcsecs were detected for the high- velocity 
components (HVC) and the low-velocity com- 
ponents (LVC), respectively. The differences 
of the derived electron density and time varia- 
tion between the two components suggest that 
these components have distinct origins, e.g., jet 
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and disc wind. The positional displacement in 
the LVC is consistent with the LVC originating 
from a disc wind. 

The forbidden lines have only blue-shifted 
components, as do those in many other T Tauri 
stars, and this is usually explained by the ob- 
scuration of the red-shifted flow by a circum- 
stellar disc. On the other hand, the Ha outflow 
has both blue-shifted and red-shifted compo- 
nents. Such a difference suggests the presence 
of a disc gap with an outer radius of 3-4 au. The 
infrared spectral energy distribution is consis- 
tent with the presence of a gap on this scale. 
This gap could be induced by a unseen com- 
panion such as a young planet. 
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Table 1: The equivalent widths of the emission lines. 
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a References - Lago & Penston (1982), Osterbrock, Tran, & Veilleux (1992) 

b Uncertainty of the measurement is 0.1 A for Ha, and 0.03-0.05 A for the others. 

c The equivalent width was not measured because of the difficulty of the removal of the strong 

Ha wing. 

d The equivalent widths of the absorption are displayed. 



Table 2: [SII] equivalent widths & ratios 



Line Component Equivalent Width (A) 

25 Aug 96 27 Jun 97 2 Jul 99 
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Note: the dividing velocity for the high velocity component (HVC) and low velocity compo- 
nent (LVC) is defined as —95 km s _1 in the stellar rest frame, as that defined by Hamann 
(1994). Uncertainty of the measured equivalent widths is about 0.02 A. 
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